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Abstract 

We investigate localization properties in the highly anisotropic and intrin- 
sically disordered layered material, which is analogous to high-T c cuprates. 
By varying the anisotropy of the system which is parameterized by the in- 
terlayer hopping tp, we find a crossover from two-dimensional (2D) to three- 
dimensional (3D) behavior at a critical hopping amplitude tp c , where a mobil- 
ity edge starts to appear. We show that below the mobility edge, anisotropic 
localization effect may exist for a finite size system, when the afr-plane local- 
ization length is longer than the system size and the c-axis localization length 
is shorter than the system size. Nevertheless, we argue that such anisotropic 
localization can not account for the "semiconductor" like behavior of the c- 

axis resistivity of high T c cuprates. 

PACS numbers: 71.30.+h, 74.25.Fy, 74.80.Dm 
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It has been controversial whether anisotropic localization could exist in intrinsically dis- 
ordered layered systems such as the high-T c cuprates There are two interesting theo- 
retical issues: first, whether it is possible at all for a three-dimensional (3D) system to be 
localized in one direction and extended in the other two directions; secondly, if it is possible, 
whether this anisotropic localization effect can account for the c-axis resistivity in high T c 
cuprates. In this paper, we address both issues. 

According to the scaling theory of metal- insulator transition ||, for an infinite 2D system, 
no matter how weak the disorder is, there is no real metallic behavior, the system is localized. 
In the presence of inelastic scattering, however, when the length scale of inelastic scattering 
L in is smaller than the localization length £, electron loses its phase coherence before it can 
finish the back scattering process, thus localization can not be realized. This is the so called 
finite size effect, the system has the same physical property as if it has a size of order of L in . 
In a physical system, the inelastic scattering length L« n is a function of temperature which 
decreases as temperature increases, thus effectively, one can mimic the effect of temperature 
by changing the cutoff length of the system. If the localization length is longer than the cutoff 
length, the system would have metallic behavior; otherwise, it would behave as an insulator. 
This is the reason why two-dimensional electron systems in semiconductor heterostructures 
have very high mobility. One expects that anisotropic transport behavior (such as in high 
T c cuprates) could arise in a finite-size (or at finite temperature) anisotropic system (which 
has different localization lengths along different directions) even if the system is below the 
mobility edge. 

While a 2D system is always localized, it is known that there is a mobility edge in 3D 
system, where metal insulator transition occurs ||. We model the highly anisotropic and 
disordered layered material by an anisotropic tight binding lattice, which is parameterized 
by hopping t (which is set to be 1 throughout this paper) in the afe-plane and tp along the 
c-axis. In the limit of tp = 0, the system is equivalent to a set of disconnected 2D planes, 
and therefore has no mobility edge. While at tp/t = 1, the model describes an isotropic 
3D system, there is a well defined mobility edge, whose position is in principle a function 
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of disorder. It is interesting to investigate the localization properties in the intermediate 
regime (0 < tp/t < 1). Previous work ||J] on this model has shown that the mobility edge 
immediately emerges once a finite tp is introduced. However, in the presence of strong 
disorders, the situation is more complicated. One expects that the mobility edge appears 
only when tp exceeds a critical value tp c which is a function of the disorder. Equivalently, if 
one fixes the energy to be at the band center E = 0, and changes tp, one expects a metal- 
insulator transition at tp = tp c , since E = is where the extended states start to form once 
the mobility edge appears. This effect is somewhat similar to the effect studied by Clark 
et al ||, in which they find a critical hopping amplitude between two coupled Luttinger 
Chains, which separates coherent from incoherent tunneling between two chains. 

In this paper, we systematically investigate the possible anisotropic localization effect in 
layered materials. The conductances for a cubic sample (along both c-axis and ao-plane) 
are calculated by using the Landauer-Biittiker formula. Our numerical results show that 
upon increasing the strength of interplane hopping through a critical tp c , c-axis conductance 
changes from semiconductor like to metallic, provided that the aft-plane is sufficiently disor- 
dered (characterized by disorder strength W), yet metallic. We attribute this phenomenon 
to the highly anisotropic localization effect - a metal-insulator transition along c-axis for a 
finite-size sample. We also calculate the finite-size localization length along the c-axis using 
the transfer matrix technique. By carrying out a finite-size scaling analysis, we found a 2D 
to 3D crossover where a mobility edge starts to emerge when c-axis hopping is above the 
same critical tp c . For highly anisotropic system (such as tp < 0.3t), the critical tp c could 
be determined from a universal relation W/tp c ~ 30 — 40. We discuss the relevance of this 
anisotropic localization effect to the anisotropic transport phenomena in high-T c cuprates. 

In the following, we briefly outline our model and technique to calculate the conductance. 
We model our anisotropic three-dimensional system in a cubic geometry with a finite width 
(M) square lattice with nearest neighbor hopping. The disorder potential is modeled by the 
on-site white-noise potential Vij m denote the column and chain indices, m denotes the 
plane index) ranging from — W/2 to W/2. The interplane hopping along c-axis tp is different 
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from the a6-plane hopping amplitude t. We calculate the conductance in both ao-plane and 
c-axis, for different degree of anisotropy levels tp/t ~ — 0.5. The Hamiltonian of this 
system can be written as: 

H = Y^ Vijm\ijm >< ijm\ (1) 

ijm 

— X] [\v m klm\ + c.c] 

m <ij;kl> 

~~ tp2_j \\^j m ^ m + 1| + c - c -] > 

ij m 

where < ij; kl > indicates nearest neighbors on the mth plane. From now on, we choose 
afe-plane hopping amplitude t as the energy unit (t = 1 throughout the paper). The dc 
conductance at Fermi energy E for a cubic sample is calculated using the Landauer-Buttiker 
formula 0: 

g = jTtT\E)T(E) } (2) 

where T(E) is the transmission matrix which is related to the retarded Greens function 
G + (0,L) across the disordered sample (along direction LinaMxMxI sample). The 
retarded single particle Greens function can be calculated using the recursive Greens- function 
algorithm [^]. 

We present in Fig.l, the calculated resistance (1,000 sample ensemble average) as function 
of the inverse of the sample size (which could mimic the effect of temperature) for both aft- 
plane and c-axis, with different strength of disorder W. All the calculations are at the band 
center, where energy E = 0. From Fig.l, we can see that for W = 2,4,6, and 8, the in- 
plane resistance increases with 1/M, displaying a metallic behavior, even more so (the slope 
increases) as the hopping amplitude tp in the c-axis increases. While the c-axis resistance 
may change from metallic to insulating by changing tp. For example, for W = 2, the c-axis 
resistance R is metallic-like for tp > 0.05, while for tp < 0.02, R increases with sample size 
M, the system appears to be localized along c direction. In this case, there seems to be a 
critical hopping amplitude tp c ~ 0.03 — 0.04 that separates the localized and the metallic 
conduction along the c-axis. As we go to higher degree of disorder, tp c increases. This 
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anisotropic localization effect (the system is only localized along c-axis) does not violate 
the scaling theory of localization || which states that there only exists a single mobility 
edge for any disordered system even if it is highly anisotropic (as shown in recent numerical 
calculations 0), i.e. for large enough system, ao-plane would undergo the metal-insulator 
transition at the same tp c as along the c-axis. Below the mobility edge, the localization 
length in the aft-plane is extremely large, thus the system could have metallic behavior 
along ao-plane if the sample size is smaller than the localization length. Nevertheless, if 
the system is big enough, we expect the a6-plane would have insulating behavior |J. It is 
worth mentioning that above the mobility edge (when tp > tp c ) , the calculated conductance 
in both a6-plane and c-axis increase linearly with the increasing sample size M, indicating 
typical 3D Drude conductance behavior. 

The finite size correction to the conductivities of an anisotropic system has been calcu- 
lated by Wolfle and Bhatt, from evaluating the maximally crossed diagrams by assum- 



ing both anisotropic single particle energy dispersions and anisotropic impurity scatterings. 
Their results show that backscattering gives rise to a diffusive correction to the metallic con- 
ductivity, and that the effects of anisotropy near the weak coupling (metallic) fixed point can 
be completely absorbed into an anisotropic diffusion coefficient. More recent work by Szott 
et al has microscopically derived the anisotropic diffusion constants from the anisotropic 
impurity scattering rates for the superlattice semiconductor, and concluded that the con- 



ductivity for superlattices still has the diffusive form. jlT] Nevertheless, these above work did 
not address the 2D-3D transition found in our numerical calculations. To be more specific, 
it is shown that the contribution from backscattering (weak localization correction) reduces 
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where / is the mean free path. Therefore, we expect that in our anisotropic layered system, 
the weak localization correction would have a crossover from logarithmic behavior to 1/M 



as tp increases from (pure 2D system) to 1 (isotropic 3D system). On the other hand, 
our numerical data show that above the mobility edge (when tp > tp c ), the calculated 
conductance in both afe-plane and c-axis increase linearly with the increasing sample size 
M, exhibiting a typical 3D Drude conductance. For large sample sizes, the breaking away 
from linear behavior in our data compare favorably to a negative weak localization correction 
in the conductivity. However, strong finite size effect and statistical fluctuations do not allow 
us to make a quantitative comparison with the analytic results of Refs |TI| and |JTT . 



In order to study the transition between 2D and 3D behaviors by tuning tp in our 
anisotropic tight-binding model, we need to calculate the localization length directly. We 
now introduce the technique to calculate the localization length in a long square bar geometry 
(M x M x L) along the c-axis. For a specific energy E, a transfer matrix T m can be easily 
set up mapping the wavefunction amplitudes at plane m — 1 and m to those at plane m + 1, 



i.e. 
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where H m is the Hamiltonian for the mth plane, J is a M 2 x M 2 unit matrix. Using a 



standard iteration algorithm ||13[| , we can calculate the Lyapunov exponents for the transfer 
matrix T m . The localization length Xm{E) for energy E at finite width M is then given 
by the inverse of the smallest Lyapunov exponent. We choose sample length L{m 1000) to 
be long enough so that the self-averaging effect automatically takes care of the ensemble 
statistical fluctuations. In our calculation, we fix the energy at E = and change the 
hopping amplitude tp for sample size M = 4 ~ 16. 

To obtain the thermodynamic properties of the system, we carry out a finite-size scaling 
analysis of the renormalized localization length Xm/M, which satisfies a simple scaling law. 
As shown in Fig. 2, all Xm/M data for different disorder and hopping amplitude tp fall on 
the same curve which can be written as: 



X M (W,tp)/M = mw,tp)/M] 
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(5) 



where tp) is a characteristic length which depends on W (or energy E) but not on 
M. Similar to the scaling analysis for an isotropic 3D system the scaling function 



has two branches. The scaling parameter associated with the lower branch corresponds to 
the thermodynamic localization length, i.e. £ = Aqq. The one associated with the upper 



branch belongs to the extended regime and identifies with the resistivity of the system |TJ 
The scaling parameter should diverge at the critical point which separates the two branches 
in the scaling function indicating a metal-insulator transition. The inset of Fig. 2 plot the 
scaling parameter £ as a function of W/tp for all the W and tp shown in the scaling curve. 
We can see that for the highly anisotropic system we considered (tp/t < 0.5), the scaling 
parameter is a universal function of W/tp. At W/tp ~ 30 — 40, £ diverges, showing a typical 
Anderson metal-insulator transition. We now compare to Fig.l, where the critical hopping 
amplitude separating the metallic and localized regimes for W = 2,4,6,8 all seem to be 
pinned at W/tp c ~ 30 — 40. Thus we argue for this highly anisotropic tight-binding model 
{tp < 0.5), a typical metal-insulator transition occurs at certain tp c (W), where W/tp c is a 
universal number independent of hopping tp and on-site disorder W. 

Here we should emphasize that the relation W/tp c ~ 30 — 40 is only true for highly 
anisotropic system where tp « t. (We fix t = 1 in this paper.) However, for systems with 
modest anisotropy, namely 0.5 < tp/t < 1, we do not expect that the relation W/tp c ~ 30 — 
40 to be relevant. In fact, we find that W/tp c ~ 16 for isotropic system (i.e. tp = t), which 



is consistent with previous works |4,14|. For the strong anisotropic system we have studied, 



where tp/t < 0.5, the system is closer to 2D, the mobility edge is mostly determined by the 
small amount of interlayer hopping tp, while the in-plane hopping t plays a much weaker 
role in determining the critical point of W/tp c other than setting the energy scale of both 
W and tp. As the system becomes less anisotropic, tp and t become compatible, therefore 
should both determine the mobility edge. Thus we find for tp/t < 0.5, the mobility edge 
appears at W/tp c ~ 30 — 40 and W/t ~ 2 — 8; when anisotropy decreases, at 0.5 < tp/t < 1, 
it is only natural that these two ratios merge closer, and eventually, at tp = t, we have 
W/tp c = W/t ~ 16. Another way to understand this is to write the critical ratio W/tp c as 



function of tp/t which characterizes the anisotropy of the system, we find W/tp c = 30 — 40 
at tp/t -C 1, and W/tp c = 16 at tp/t = 1. We expect a smooth crossover in the intermediate 
anisotropy regime of 0.5 < tp/t < 1. Our results show that W/tp c is overall t-dependent, 
while in the strong anisotropic limit of tp/t < 0.5, the t-dependence of W/tp c becomes much 
weaker. 

In conclusion, we observe a crossover from 2D to 3D behavior in the highly anisotropic 
layered system when the interlayer hopping amplitude exceeds a critical value. This crossover 
behavior of localization from 2D to 3D has been studied by various researchers, such as in Ref . 
4. However, we have investigated the non-trivial consequences of this crossover behavior by 
directly calculating the conductance of finite-size samples. Anisotropic transport properties 
(such as temperature dependence) could arise from the anisotropic localization properties. 
We have shown that even when the system is below the mobility edge, in the presence of 
inelastic scattering which provides a cutoff length scale, it is possible to observe localization 
in the c-direction (if the localization length along c-axis is shorter than the cut-off length) 
and metallic behavior in the a6-plane (if the localization length in the afr-pane is longer than 
the cut-off length) for highly anisotropic systems. 

A natural applicable system where the above anisotropic localization may occur is 
the layered high T c cuprates, where the c-axis resistivity indeed shows a "semiconductor- 
like" upturn at low temperatures for Bi 2 Sr 2 CaCu20 8 and underdoped La 2 _ x Sr x Cu0 4 
and YBa 2 Cu 3 6+x , while the upturns disappear upon doping. Nevertheless, optical far- 
infrared and Raman experiment [[HJ suggest that the c-axis scattering rate l/r c for typical 



YBa 2 Cu306+x is of the order of lOOmeV, while the c-axis hopping amplitude tp ~ lOmeV. 
Thus it appears that the ratio of disorder (should not be larger than l/r c if there is any) and 
c-axis hopping W/tp < W/tp c . Moreover, the typical afe-plane hopping energy is of the order 
of lOOmeV, leaving W/t < 1. Therefore, these cuprates are well above the mobility edge, 
and the c-axis transport should be on the metallic side. Details of the c-axis scattering rate 
and the hopping rate for each individual cuprate at different dopings are listed in Ref. [TIB 



Furthermore, unlike the simple model we present here, there is a strong afe-plane fluctuation 



which increases with temperature, unless at very low T, this planar fluctuation will destroy 
the phase coherence so that back scattering can not occur 0. Therefore, we argue that the 
c-axis resistivity of high-T c cuprates can not be associated with the anisotropic localization 
phenomena we discussed above. An reasonable understanding of the c-axis resistivity data 
is given in the recent work of Ref. |TB[ . Nevertheless, this does not exclude the possibility 
that localization (either c-axis or a6-plane) may occur in some extremely dirty samples or 
films or at doping concentration in closely proximity of the insulating phase of cuprates. 

We acknowledge stimulating conversations with A.J. Leggett, P. Philips, K. Levin, J. 
Miller, and Q.M. Li. This work is supported by NSF-DMR-91-20000 and 94-16926 through 
the Science and Technology Center of Superconductivity. 
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FIGURES 

FIG. 1. Sample size dependent resistance for a cubic sample for different disorder strength (a) 
W = 2; (b) W = 4; (c) W = 6; (d) W = 8. 

FIG. 2. The scaling function, Xm/M vs £/M. The £(W,tp) shown in the inset turns out to be 
a universal function of W/ tp. 
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